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We propose a model for the non-Fermi behavior in the proximity of the quantum phase transition 
induced by the strong polarization of the electrons due to local magnetic moments. The self- 
consistent Renormalization Group methods have been used to calculate the temperature dependence 
of the electrical resistivity and specific heat. The y 5 / 3 dependence of resistivity and the TlnT 
dependence of the specific heat show that the magnetic impurities drive a ferromagnetic quantum 
phase transition and near the critical point the system present a non-Fermi liquid behavior. The 
model is in good agreement with the experimental data obtained for Ni^Pdi-x alloy. 



I. INTRODUCTION 

Deviations from the standard Fermi liquid descrip- 
tion in heavy fermion systems (HFS) have been asso- 
ciated to the proximity of the quantum phase transi- 
tion (QPT). Denned as phase transitions at T = 0, 
QPT's are usually driven by quantum fluctuations con- 
trolled by a non-thermal parameter, namely by impu- 
rities, pressure or magnetic fields. Most of these low 
temperature magnetic QPT's are from paramagnetic to 
antiferromagnetic state, but recently HFS undergoing a 
paramagnetic to ferromagnetic phase transitions were 
identified. Example of HFS which are close to fer- 
romagnetic order are Tb.i_ x U x Cu2Si2,|il MnSi,[j| and 
Ni^Pdi-a;. 4] The itinerant-electron ferromagnetic state 
in these materials is induced by pressure (MnSi) or im- 
purities (Th 1 _ a .U 2; Cu2Si2 and Ni a Pdi_ x ). 

Thi_ :c U a ;Cu2Si2 for concentrations x > 0.15 presents 
ferromagnetic order. Experimental data in these com- 
pounds show a non- Fermi liquid behavior for C /T oc In T 
and a critical temperature, T c , as low as 12K for the 
x = 0.15 sample. However, the magnetic susceptibil- 
ity data present a field dependance even at temperatures 
bellow T c , a behavior characteristic for weak itinerant 
ferromagnets. 

In the case of MnSi experimental results showed the ex- 
istence of a first-order phase transition induced by pres- 
sure, the critical temperature decreasing towards abso- 
lute zero at a pressure value p c ~ 14.6 kbar. For p > p c 
a non- Fermi liquid behavior of the resistivity as function 
of temperature is reported, p ~ T.|3( The phenomenolog- 
ical theory Q based on the interaction of the electrons 
with the overdamped spin fluctuations at low energies, 
can explain the occurrence of the non-Fermi state. The 
quantum transport anomalies of the itinerant-electron 
ferromagnetic state have been discussed also by Bclitz, 
Kirkpatrik, Narayanan and Vojta (BKNV) the non- 
Fermi behavior of the system being proved based on the 
scaling approach. The model from Ref. |{| predicts a first 
order phase transition and it is in an excellent agreement 
with the QPT driven by pressure in MnSi 

NLcPdi-z at a Ni concentration x c = 0.025 presents 
ferromagnetic order assumed to occur in the itinerant 
electron system of Pd atoms due to a strong polarization 



of these electrons by magnetic impurities (in this case 
Ni). The transition critical temperature depends on the 
Ni concentration as (x — x c ) 3 ' 4 in the critical region. In 
the same region of the phase diagram, experimental data 
reveal a T 5 / 3 dependance of the relative resistivity, Ap, 
and a TlnT dependance of the specific heat, C(T).[4( All 
these experimental data clearly identify the non-Fermi 
character of the electronic excitations close to the quan- 
tum critical point (QCP). These results were successfully 
described in terms of phenomenological spin-fluctuation 
models, p 

In this work we extend the Doniah-Wohlfarth model, 
proposed for the explanation of the itinerant electron 
ferromagnetic state driven by impurities in Ni^Pdi-^ 
compounds, to the critical region and we calculate using 
renormalization group methods the temperature depen- 
dence of the resistivity and the specific heat. We will 
apply the Hertz-Millis H @ version of the Renormal- 
ization Group method (RNG) to the Doniah-Wohlfarth 
model taking into consideration the quantum effects at 
finite temperature and extract the effects of the spin fluc- 
tuations on the system. The self-consistent renormaliza- 
tion group theory, given by Moriya ^} will be used to 
calculate the temperature dependence of the electrical 
resistivity. 



II. MODEL 

We consider that the fluctuations in the magnetization 
in the critical region are given by the action: 

s eff = si% + si% , (1) 

where 

s Sf = rt) w«)i a ( 2 ) 
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and 
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Here we introduce the notation q = (q, ui n ), u> n being the 
bosonic Matsubara frequencies, and 



E = k * T E 



d d q 

(2n) d 



In Eqs. ©-Q x(q) 1S the fluctuation propagator and u 
is the coupling constant. 

In the following we consider that the spin impurity has 
a very strong polarization effect on the electrons and at a 
critical concentration x — x c a new phase, which is in fact 
a ferromagnetic phase, can be reached. This model has 
been proposed first by Donniach and Wolhfarth [?J using 
a single impurity approximation. The susceptibility of 
the polarized electrons was given (see Ref. Q) as 



xo(q,^) 



1 



1 + 2 



J 2 R> 
JR-u 



xo(q,w) 



(4) 



where Xo(q, is the susceptibility of the electronic sys- 
tem. In Eq. Q J is the exchange interaction between 
electrons and localized spins and / is the interaction be- 
tween electrons. The parameters R and B! have been 
calculated as 



R = N £ ^ ~ 



rikt) 



and 



R' 



-x (S z 



(5) 
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x being the impurity concentration of the magnetic mo- 
ments with spin S. 

The dynamic susceptibility xo(q, has the form ^jj] 



xo x (q,w) -x _1 (o,o) 



1 - Dq 2 + iC- 



(7) 



where D and C are constants. We approximate x(q, ui) 
from Eq. Q as 



x(q,w) 



i 



x -*(q )W )-J-2^ 
a results which based on Eq. Q can be written as 

x(q,w) = 



5 (x) + aq 2 - if^ 



(8) 



(9) 



where a and Y are constants, and 5o(x) is given by 

*o(x) = Xo 1 I 2^ . (10) 

The parameter Sq(x), with a linear dependance on the 
impurity concentration, measures the distance from the 
QCP. 

The model is valid only for systems in which the lo- 
cal moments give a strongly polarization of the itinerant 



electrons and this is typically for metals as Pd which is 
paramagnetic and the Stoner criterium cannot be sat- 
isfied only due to the electron-electron interaction. For 
this particular system we do not expect localization ef- 
fects, excepting the case of nonmagnetic impurities, but 
this cannot drive the system in the ordered phase. 

This model, described by Eqs. QJ, © and ©, can be 
treated using the Renormalization Group method, in the 
versionproposed by Hertz and Millis 0- Following 
Refs. |3 01 we perform the standard scaling k — ► k'/b, 
uj n — > tu' n /b z , z being the dynamic critical exponent. We 
obtain the following flow equations: [TT| 



dT(l) 
dl 



dY(l) 



dl 



zT(l) , 



(3-z)r(0, 



dS(l) 
dl 



28(1) + 2(n + 2)/m(Z) , 



(11) 



(12) 



(13) 



du(l) 
dl 



[A-(d+z)]u(l) - (n + 8)f 2 u 2 (l) , (14) 



where d is the spatial dimension, h are functions char- 
acteristic for the model (see Ref. Q) and I = In b is the 
scaling variable. Additionally, the system free energy will 
scale as: 



dF(l) 
dl 



(d + z)F(l) + h 



(15) 



/3 being again a characteristic function of the model. 



III. SPECIFIC HEAT AND RESISTIVITY 

The evaluation of the renormalized free energy F[T(l)] 
based on Eq. (|15fl will lead to the temperature depen- 
dence of the specific heat which by definition can be cal- 
culated as the second derivative of the free energy with 
respect to the temperature, C(T) = -T (d 2 F / dT 2 ) . In 
general there will be two distinct contributions to the 
renormalized free energy, associated to a quantum do- 
main, T(l) <C 1, and to a classical one, T(Z) 3> 1. For 
more details on the calculation of the free energy see Refs. 

However, as the form of the system susceptibility 
match the corresponding form for a ferromagnetic sys- 
tem, if we consider the d = 3 case, the specific heat is 
obtained as 



C(T) = l0 T + 7l T InT, 



(16) 



7o and 71 being constants, a result which clearly show 
that the behavior of the considered system is non-Fermi, 
as corrections to the linear temperature dependance of 
the specific heat are logarithmic. This result is in agree- 
ment with the experimental data presented in Ref. Q. 
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We have to mention that a similar behavior was obtained 
using RG method for a system in the proximity of the 
Lifshitz quantum critical poin t. [T3| Recent experimental 
data in silicon MOSFETs[13| suggested a QPT to a fer- 
romagnetically ordered state in d — 2. An analysis of 
the s peci fic heat behavior for the d — 2 case was done in 
Ref. [nj suggesting a different temperature dependance 
of the specific heat. 



In order to calculate the temperature dependence of 
the resistivity we apply the self-consistent theory of fluc- 
tuations to the action given by Eq. Q . This can be done 
in the version of 1/n expansion (n being the number of 
components of the bosonic field 0) applied to the (p 4 ac- 
tion. Using the approximation 4 | ~ 2 < \<fi 2 \ > \<fi 2 \ the 
renormalized parameter 6(x) can be calculated from the 
following self-consistent equation: 



(2tt) 3 8{x) + k 2 + 



(17) 



The summation over the bosonic Matsubare frequencies on the second term in the right hand side (rhs) of Eq. i|17|) 
can be performed analytically leading to the following expression: 



S(x) = 5, 



d 3 k 

(2tt) 3 
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2k B Tj (S + h 2 ) 2 + (^y 
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The temperature dependence of the QPT parameter 5(x) can be extracted if we consider Eq. (|18J) at the QCP in 
order to eliminate the bare QPT parameter Sq(x). Accordingly, we obtain 



S{x) = -(- + 1 



2 V2 



d 3 k 

(2tt) 3 

ri 3 k 

(^o 3 



rfc 



duj 



coth 



rfc 



2k B T 



rfc 



rfc 



{d{x) + k 2 ) 2 + (^y 

Tk 



(S(x) + k 2 ) 2 + (^) 2 fc 4 + (^) 2 



(19) 



For T = we will show that 5 < T holds. In this 
approximation the first integral in the rhs of Eq. i|19|) 
becomes: 



h = 2k B T 



d A k 



k B T i uj 

duo n 



(2tt) 3 Jo u A 2 + (fj:) 



2 ' 
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where A = 8{x) + k 2 . Performing the integral over u> we 
get: 

4/3 
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where y = k B T/Tk 3 ; y c , associated with the upper crit- 
ical wave- vector k c , has been substituted by infinity. C\ 
is a constant. The second integral from the rhs of Eq. 
02), I 2 , is: 



rfc 



do; 
2^ 



rfc 



rfc 



^ 2 + (rt) 



fc4 + (r%) 



(22) 



and can be performed if one introduce a new variable 
y = k/8 1 / 2 with an upper cut-off y c — fc c [(5(a:)] _1 ' /2 . The 
final result can be express as: 



-C 2 5{x) 



(23) 



Based on Eqs. J2U an d 1221 w ^ obtain for 8{x) the 

following temperature dependance: 

S(x,T) = C(u,n) (fc s T) 4 / 3 , (24) 

which is finite for n — > oo. This result can be inverted 
in order to extract the concentration dependence of the 
critical temperature. If one consider that the tempera- 
ture can be approximated by the critical value (T c ), and 
that in the first order S(x) ~ Sq(x) we have 



T c (x) ~ (x - x c f' 4 



(25) 



a result which is in agreement with the one discussed in 
Ref. g. 

The temperature dependance of the resistivity can be 
extracted as the imaginary part of the self-energy, ob- 
tained as a result of electrons interacting with the ferro- 
magnetic fluctuations. In the one-loop approximation we 
have: 



S(k, iu> r , 



g 2 k B T Y, G ( k + q; ioj n 



iui) D(q,iuji) , 
(26) 



where g 2 is the coupling constant, G(k, iu) n ) is the elec- 
tronic Green function. 



G(k, iw n ) 



1 



iuj r . 



e(k) 



(27) 
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with u> n = (2n + l)fcsT, and Z?(q, iu>i) has the same form 
with x(q, 0) given by Eq. whit 5 (x) replaced by S(x). 
Performing the summation over lo\ in Eq. (j26|) we obtain: 



E(k,iw n ) = -ng 2 V / — 

„ J -To 11 



dz Yq 



q j-r 9 A 2 + ^ 
x [n B (z) + n F {z + uj)} 5[lo + z - e(k + q)] , (28) 

where ub(z) is the Bose function and np(z) is the Fermi 
function. For \z\ <Twe approximate n_e(z) + np(z) ~ 
T/z and performing the analytical continuation iia n — > 
lu + irj we calculate 



-7m S fl (/c F ,0) 



7rg 



E 



<J[z - e(k + q)] (29) 



which leads to the following expression for the imaginary 
part of the self-energy 



Im T, R (k F ,0) = const 



{k B Tf 
5{T) 



(30) 



Using now Eq. I|24|l we obtain the temperature depen- 
dence of the scattering time l/r e ff = —IrriE R as: 



r eff 



const (ksT) 



5/3 



(31) 



which gives for the temperature dependent resistivity 
Ap(T) = p(T) - p(0) the following behavior 

Ap(T) ~ (fc s T) 5 / 3 , (32) 

in agreement with the experimental data from Ref. 

IV. CONCLUSION 

In conclusion, we presented a model for the ferromag- 
netic QPT driven by the magnetic impurities which po- 
larize the Fermi liquid close to the Stoner instability. The 



result can be regarded as a generalization of the Doniach 
and Wohlfarth mean-field model for the case of QPT. 
As a general result we conclude that the system presents 
a non- Fermi behavior in the critical region around a QCP. 
The non- Fermi character of the system is sustained by a 
T In T behavior of the specific heat correction term, and 
by an electrical resistivity which presents a T 5 / 3 depen- 
dence. Both these results are in good agreement with the 
experimental data reported in Ref. [4J. A similar tem- 
perature dependance of the resistivity was obtained by 
Mathon[l4( using a simple molecular field theory. How- 
ever, despite a good agreement between the calculated 
and experimental values of the resistivity, Mathon calcu- 
lations are not able to explain the non-Fermi behavior of 
the specific heat and do not take into account the quan- 
tum effects observed in Ni^Pdi-^, making the model in- 
adequate for the proper description of the QPT in this 
particular system. A phenomenological description of 
the QPT was also made by Lonzarich. (6j We also explain 
the concentration dependance of the critical temperature, 
our result being in good agreement with the experimental 
data. 



We mention that this model is different from the 
BKNV p| model, in which the spin susceptibility for the 
d — 3 case is considered as x(g, lu = 0) = Xo+q 2 ^(pf/q)- 
This form of xio) h as been carefully analyzed by Mil- 
lis 15] and at the present time new experimental data 
are needed for a confirmation of this spin susceptibility. 
Recently, Belitz and Kirkpatrick 0] have been reconsid- 
ered the QPT in the clean itinerant-electron ferromagnet. 
The coupling of the order parameter fluctuations to the 
soft fermionic fluctuations lead to a theory which is very 
different than theories based on the Hertz-Millis model. 
The main point in the new version of the BKNV theory 
is that the fluctuations can change the first order phase 
transition in a second order one. However, the occurrence 
of two dynamical exponents z and z for the two kind of 
fluctuations makes the two theories very different, even 
if in any case the mean field behavior can explain the 
experimental data. 
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